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In this contribution we present results from quenched QCD simulations with the parameterized fixed-point 
(FP) and the chirally improved (CI) Dirac operator. Both these operators are approximate solutions of the 
Ginsparg- Wilson equation and have good chiral properties. We focus our discussion on observables sensitive to 
chirality. In particular we explore pion masses down to 210 MeV in light hadron spectroscopy, quenched chiral 
logs, the pion decay constant and the pion scattering length. We discuss finite volume effects, scaling properties 
of the FP and CI operators and performance issues in their numerical implementation. 



1. Introductory remarks 

Fixed point and chirally improved fer- 

mions provide an approach to chiral symme- 
try on the lattice based on the Ginsparg- Wilson 
equation 0. In an actual implementation both 
the FP and CI fermions give rise to an approx- 
imate solution of this equation and it has to be 
tested how well a chiral fermion is described by 
such an approximation. Here we report on our re- 
sults from quenched QCD calculations focussing 
in particular on observables sensitive to chiral 
symmetry (for a recent general review of results 
with chiral actions see ||]). 

In another contribution to these proceedings 
|9j we gave a more general overview of our cal- 
culations and also a brief introduction to the 
construction of the FP and GI operators {Dpp, 
Dci)- Here we would like to deepen the discus- 
sion and present selected topics in more detail. 
These topics include light hadron spectroscopy. 



*This contribution is based on parallel talks and posters 
presented by C. Gattringer, S. Hauswirth, K. Holland, 
K.J. Juge, C.B. Lang and S. Schaefer at LATTICE 2002. 



the quenched chiral logarithm, the pion decay 
constant and a preliminary study of the 1 = 2 
pion scattering length. 

Both the FP and the CI operators make use 
of the full Clifford algebra and a large number 
of gauge paths to approximate a solution of the 
Ginsparg- Wilson equation. They differ in the 
method for determining the coefficients in front 
of the individual terms of the Dirac operator. 
The FP operator is constructed from classical 
equations which determine the fixed point of a 
renormalization group transformation in QCD. 
For the CI operator the Ginsparg- Wilson equa- 
tion is mapped onto a system of coupled equations 
for the coefficients of the individual terms of the 
Dirac operator and this system is then solved nu- 
merically. In our implementation of the two op- 
erators we restrict ourselves to terms essentially 
on the hypercube only. For a more detailed dis- 
cussion of the operators see [^,^,^,|[^,^. 

The parameters of our quenched QCD calcu- 
lations were chosen similar for the two operators 
such that a direct comparison of the results is pos- 
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mps/mv 


FP 


16^ X 32 


0.16 fm 


200 


0.28-0.88 


FP 


12^ X 24 


0.16 fm 


200 


0.3-0.88 


ov3 


12^ X 24 


0.16 fm 


100 


0.21-0.88 


FP 


8^ X 24 


0.16 fm 


200 


0.3-0.88 


FP 


12^ X 24 


0.10 fm 


200 


0.34-0.89 


FP 


16^ X 32 


0.08 fm 


100 


0.36-0.89 


CI 


16^ X 32 


0.15 fm 


100 


0.38-0.85 


CI 


12^ X 24 


0.15 fm 


100 


0.36-0.85 


CI 


8^ X 24 


0.15 fm 


200 


0.33-0.85 


CI 


16^ X 32 


0.10 fm 


100 


0.33-0.92 


CI 


12^ X 24 


0.10 fm 


100 


0.32-0.92 


CI 


16^ X 32 


0.08 fm 


100 


0.40-0.95 



Table 1 

Statistics for the FP and CI Dirac operator to- 
gether with the lattice spacing measured with the 
Sommer parameter and the range of pseudoscalar 
to vector mass ratios we worked at. 



sible. For both operators we used lattices of size 
8^ X 24, 12^ X 24, 16^ x 32 with lattice spacings 
of approximately a — 0.08 fm, 0.10 fm and 0.15 
fm. Thus, we were able to study finite volume ef- 
fects (constant a) as well as discretization effects 
(constant physical volume) . The FP operator was 
used together with quenched configurations from 
a parameterized perfect gauge action |]lOj subse- 
quently smoothened with renormalization group 
inspired smearing The latter can be viewed 
as part of the parameterization of the FP opera- 
tor. We also employed a version of the FP oper- 
ator (ov3) augmented with overlap steps It 
is constructed from the parameterized FP opera- 
tor with three terms of the Legendre expansion of 
the overlap. For the chirally improved operator 
we generated the quenched gauge configurations 
with the Liischer-Weisz gauge action |12 and ap- 
plied one step of HYP blocking |l^ to improve the 
properties of the Dirac operator. An overview of 
the statistics is given in Table 0. There we also 
give the results for the lattice spacing as deter- 
mined from the Sommer parameter [ p^Jl4| and 
the range of pseudoscalar to vector meson mass 
ratios we have worked at. The smallest pion mass 



we have worked at is 210 MeV for the FP opera- 
tor (160 MeV with ov3) and 240 MeV for the CI 
operator. We did not encounter problems with 
exceptional configurations for the listed range of 
pion masses and we expect that we can reach pion 
masses of 200 MeV for both operators without 
having to go to exceedingly small lattice spacings. 

Our hadron propagators were created by 
smeared sources. For the FP operator we fixed 
the gauge configurations to Coulomb gauge and 
used Gaussian sources. For the CI operator we 
employed Jacobi smearing JTsI and no gauge fix- 
ing was necessary. 

At every fermion offset of the hypercube a large 
number of gauge paths contribute in our Dirac 
operators. By precalculating and storing some 
basic gauge link products the construction of the 
corresponding sparse matrix could be speeded up 
significantly. Every offset on the hypercube is 
represented by a color-Dirac matrix whose size 
is 12 X 12 independently of the number of gauge 
paths to this offset. When comparing the cost 
to e.g. the Wilson operator (without clover term) 
the number of fioating point operations in the ba- 
sic matrix-vector multiplication is increased by a 
factor of w 36. However, since the bottleneck is 
the access to memory and we do more operations 
with each SU(3) matrix we fetch from memory 
(all 16 entries in the Dirac space are used, Wil- 
son uses only 4) the actual number might be quite 
different. Furthermore we observed Q that some 
numerical operations such as the computation of 
low lying eigenvalues converge considerably faster 
for our Dirac operators than for Wilson's opera- 
tor. 

Our calculations were done on the Hitachi 
SR8000 at the Leibniz Rechenzentrum in Mu- 
nich and the production runs were typically per- 
formed on 8 to 16 nodes each consisting of 8 CPUs 
with shared memory. The Dirac operator was dis- 
tributed among the nodes and the parallelization 
was done using MPI. 

2. Pion masses 

2.1. Different correlators for the pion 

In order to extract the pion we experimented 
with different operators which have overlap with 
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this channel. In Fig. |l] we compare the effec- 
tive mass plots at different bare quark masses ob- 
tained for three combinations of sources and sinks 
constructed from the point-like pseudoscalar and 
axial- vector densities. The mass plateaus are 
quite well pronounced although for the third cor- 
relator the signal is noisier as it is expected for a 
correlator with additional derivatives. 

The masses corresponding to the 2-point func- 
tions were extracted from overall fits to cosh- 
behaviour covering lattice distances between 
(^min, imax)- Thcsc parameters were fixed by in- 
specting the effective mass plots and checking 
the of the fit. In Fig. ^ we compare the 
masses from the three correlators and find that 
they agree within errors. 

2.2. Zero mode effects in pion correlators 

In the chiral region the pion correlators suf- 
fer from a topological finite size effect specific to 
the quenched approximation [ p^[[r^ . Zero modes 
of the Dirac operator (suppressed by the fermion 
determinant in a dynamical calculation) lead to 
an unphysical increase in the pion mass at small 
quark masses. Since the abundance of zero modes 
scales only as \fV with the volume, the net effect 
goes as For small volumes, however, this 

effect is quite important. 

Removing simply the contribution of the zero 
modes from the quark propagator is a dangerous 
procedure since it might change the pion propa- 
gator in a non-local manner. Another possibility 
is to consider the difference of the pseu- 

doscalar and scalar 2-point functions 
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This method is based on the fact that for ex- 
actly chirally symmetric actions the scalar prop- 
agator has the same topological finite size effect 
as the pseudoscalar propagator, while for small 
quark masses the lightest particle in this channel 
is much heavier than the pion. This observation 
suggests the following strategy. At small quark 
masses, where Gpp is strongly distorted, we de- 
termine the pion mass from Gpp-ss- Going to- 
wards heavier quark masses the significance of the 
zero modes is suppressed and we expect a window 
where Gpp and Gpp^ss l^ad to consistent mass 
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Figure 1. Effective mass plots for three differ- 
ent correlation functions (top to bottom: {P P), 
{At At), {PdtAt )) for Dci, 16^ x 32, Plw = 8.7, 
a = 0.08 fm; they demonstrate the pion content 
of the involved operators. As expected the signal 
becomes noisier from top to bottom. 
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Figure 2. Pion masses obtained from different 



correlators {Dci, Id x 32, Plw 
0.08 fm). 
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fits. In this window and beyond it we use the 
Gpp correlator. At heavy quark masses the mass 
difference between the pseudoscalar and scalar 
masses is too small to be resolved and the fit- 
ted mass from Gpp-ss would be larger than the 
true pseudoscalar mass. Fig. || illustrates these 
expectations on a small lattice for the FP-ov3 op- 
erator, fn the top plot the circle indicates the 
window where the zero mode effects are already 
negligible, but the scalaris much heavier than the 
pion and so the correct pion mass is easily seen in 
the Gpp-ss correlator. The bottom figure illus- 
trates the danger of separating modes from the 
quark propagator itself: although at small quark 
masses one obtains the same pion mass as from 
Gpp-ss, at intermediate and large quark masses 
the pion mass obtained is incorrect. No window 
exists in this case. 

Once the lattice volume is sufficiently large the 
topological finite size effect is no longer visible. 
In Fig. ^ we show the pion mass for a box with 
spatial extent of i = 2.6 fm (16'^ x 32 lattice, 
a = 0.16 fm, FP operator). The data is fit- 
ted with a quadratic polynomial (full curve) and 
also a curve including the quenched chiral log jl^ 
(dashed curve). Although the difference between 
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Figure 3. Effect of the subtraction of the zero 
modes on the pion mass (FP-ov3 operator 6"^ x 16, 
a — 0.16 fm). The mass from Gpp is represented 
by filled circles, while we use open circles to rep- 
resent the results from the subtracted correlator 
of Eq. (|l|) (top plot), and from the explicit re- 
moval of the zero modes from the quark propaga- 
tors (bottom plot). 

these fits is not seen on the scale of the figure, 
the quality of the QxPT fit is significantly better. 
Both curves extrapolate reasonably well to zero 
with the quark mass indicating that the topo- 
logical finite size effect is negligible in this large 
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Figure 4. Pion mass (circles) and AWI mass 
(crosses) as a function of the bare quark mass. 
16^ X 32 lattice, a = 0.16 fm, FP operator. 

box. In principle one can determine the quenched 
chiral log parameter directly from the pion with 
degenerate quark masses (see e.g. for recent 
such determinations with the overlap operator) 
but in the next section we study the quenched chi- 
ral logs using a more effective method. Fig. ^ also 
includes data for the axial Ward identity (AWI) 
mass (crosses) which we will discuss later. 

2.3. Quenched chiral logs 

Let us now discuss the determination of the 
quenched chiral log parameter d with a different 
method The idea is to use not only data 

for the pion with degenerate quark masses but to 
also take into account data computed with non- 
degenerate quark masses. The prediction from 
quenched chiral perturbation theory for the de- 
pendence of the pseudoscalar meson mass mps,i2 
on the quark masses mi and m2 reads |L8| 
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Figure 5. x — y plot for the FP operator (16'^ x 32 
lattice, a = 0.16 fm). 

(2) 

with A, B, 6 and ax a priori unknown constants. 
The arbitrary scale A^ is of the order 1 GeV. The 
dependence on the constants A, B and A^ can be 
removed by forming the following cross ratio y 

Hi 



y 



2mi "ips,i2 2m2 "^ps,i2 



mi -|- m2 mpg mi + m2 m^ 



(3) 



PS, 22 



For small 5, ax and small quark masses y is ex- 
pected to behave like 



y^l + 5x + axz -t- 0{77i,S'^) 



with 
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1 / 2mim2 
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, m2 

in 

mi 



(4) 



(5) 



-mi — m2 I , (6) 



which means that S can be extracted for small 
quark masses from the slope of y as a function 
of X. Since in this analysis the quark mass en- 
ters only in ratios, for mi and m2 any definition 
can be used where the quark mass has no addi- 
tive renormalization. (The multiplicative quark 
mass renormalization cancels.) The mass 77iawi 
defined by the axial Ward identity has this prop- 
erty and was used when comparing Eq. (H) with 
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Figure 6. x — y plot for the CI operator (16'^ x 32 
lattice, a = 0.15 fm). 

the data.0 The plots are shown in Fig. ^ for the 
FP operator (16'^ x 32 lattice, a = 0.16 fm) and 
in Fig. ^ for the CI operator (16'^ x 32 lattice, 
a = 0.15 fm). Suppressing axz and other higher 
order terms in Eq. (||) we find S — 0.17(2) for the 
FP operator and S = 0.18(3) for the CI operator. 

We have not yet studied the systematical errors 
of these predictions due to cut-off or chiral sym- 
metry breaking effects (our Dirac operators sat- 
isfy the Ginsparg- Wilson relation approximately 
only). Further, the data with large negative x 
come from large ratios of the two quark masses 
which implies in our case that one of the quark 
masses is large, stretching the validity of Eq. (|^). 

3. Pion decay constant 

3.1. The residual quark mass 

For a chirally symmetric theory one may de- 
fine conserved, covariant]^ currents which do not 
require renormalization and covariant scalar and 
pseudoscalar densities whose renormalization fac- 
tors satisfy Zs ^ Zp ^ 1/Z„, Here Z„ 
is the renormalization factor of the quark mass 
which is the coefficient of the covariant scalar 
density in the Dirac operator. Although the 



^The bare quark mass in our Dirac operators lias a small, 
but non-zero additive renormalization, see Sect. 3.1. 
^with respect to chiral transformations 



overhead using such currents and densities is ex- 
pected to be relatively small we used point- 
like, non-conserved and non-covariant operators 
in this work. (The mass operator in the action is, 
however, covariant.) 

In Fig. 1^ we show the ratio 



Pit) 



{p-iO)dtA+{p^O,t)) 
(F-(0)P+(p = 0,t)) 



(7) 



where P and At are the pseudoscalar density and 
the time component of the axial-vector current, 
respectively and = ± iP^, where 1,2 are 
flavor indices. This quantity exhibits excellent 
plateaus for separations above 3 lattice spacings. 
We denote the plateau values by po- The sources 
P(0) are smeared, the sink operators unsmeared. 

For renormalized fields the ratio is just twice 
the renormalized quark mass (see also the discus- 
sion in the next section); in fact, po/2 may be 
employed to define a bare mass parameter. This 
mass has no additive renormalization: po/2 = if 
rriTr = 0. We have referred to this mass before as 
the Ward identity mass niAwi- We identify the 
residual quark mass for our Dirac operators by 
the value of the bare quark mass in the action for 
which po/2 — 0. Fig. || gives an example of our 
results, both for Dqi and Dpp (at slightly differ- 
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Figure 7. Results for the ratio p{t) for Dqi, 16^ x 
32, Plw = 8.7, a = 0.08 fm. 
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Figure 8. The plateau values for the ratio 
Po(m) = 2mAwi(fTi) vs. the bare quark masses 
for both Dirac operators Dci {Plw = 7.9, a = 
0.15 fm) and Dpp (Pfp = 3.0 , a = 0.16 fm) for 
lattice size 16'^ x 32. The linear fit indicates a 
very small residual quark mass. 



ent lattice spacings) and we find linear behaviour 
for small bare quark masses. A linear fit to these 
data and data for other parameters demonstrates 
a very small residual quark mass shift ranging for 
Dci from 0.002(1) at a = 0.15 fm to 0.000(1) 
at a = 0.08 fm, and for Dpp from -0.0006(4) at 
a = 0.16 fm to -0.0194(2) at a = 0.08 fm. The 
largest value for Dpp is likely related to the fact 
that Dpp was optimized at a = 0.16 fm and then 
also used at the smaller lattice spacing without 
readjustment. 

3.2. Za and 

Exact chiral symmetry on the lattice implies 
that the covariant conserved axial current and the 
covariant pseudoscalar density satisfy the equa- 
tion 

d^A^{x)^2mP^ix) , (8) 

which can be used in on-mass-shell Green's func- 
tions, or put differently, Eq. (||) is valid in general 
bare Green's functions up to contact terms. In 
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N';^ X Nt a = 


0.15 fm 


0.10 fm 


0.08 fm 


8^ X 24 


0.94(2) 






12^ X 24 


1.00(2) 


0.96(1) 




16^ X 32 


1.00(2) 


0.97(1) 


0.96(1) 



Table 2 

Values for the renormalization factor Za of the 
point- like axial current for Dci- 



this equation m is the bare quark mass multiply- 
ing the covariant scalar density in the Dirac oper- 
ator, Afj^ is the covariant conserved current and P 
is the covariant density. Eq. (^ is a consequence 
of bare Ward identities which follow directly from 
the path integral formulation of lattice QCD with 
Ginsparg- Wilson fermions. 

Considering Eq. (^ with such covariant and 
conserved operators we get 2m on the left hand 
side. Replacing the conserved covariant current 
on the r.h.s. by ZAAt, where At is the naive 
point-like current, we can calculate Za in terms 
of the correlators and m. It can be shown jj] 
that, for 2R = 1 in the Ginsparg- Wilson rela- 
tion as it is the case for Dci, in these correla- 
tors the covariant pseudoscalar density can be re- 
placed by the naive point-like density divided by 
(1 - m/2)(l - m^/4). This factor goes to 1 in 
the chiral or in the continuum limit. At the end 
the renormalization factor Za of the point-like 
current can be obtained by measuring correlators 
of the point-like operators. Up to factors (~ 1) 
mentioned above Za = 2m/ po, where po is the 
plateau of the ratio in Eq. (^ discussed in the 
previous section. We collected the corresponding 
numbers in Table |^. These numbers might be 
distorted by the small chiral symmetry breaking 
effects in Dci- The trick of replacing the covari- 
ant densities by their naive point-like versions is 
modified for the case when 2R is a non-trivial lo- 
cal operator [|j as it is the case for Dpp- For the 
FP action the corresponding correlators are not 
yet measured. 

Similar considerations can be used to con- 
nect with the correlator of point-like, naive 
pseudoscalar densities. Start with the correlator 
^-(P(a?, i)P(O)) where P is the covariant den- 
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Figure 9. Results for U for Dci, 16^ x 32, 
Plw = 8.7, a = 0.08 fm together with fits hn- 
ear or quadratic in a m. 



Figure 10. The (quadratic) extrapolations of 
to the chiral limit for different lattice spacings ex- 
hibit simple scaling behaviour (results for Dqi)- 



sity entering the AWL Saturating this correlator 
by the pion intermediate state for large time sep- 
aration, using Eq. (H) and the definition of the 
pion decay constant^ one obtains 

~ 8m'^ \ ) ■ 

As before, the covariant pseudoscalar density in 
Eq. ^ can be replaced by the point-like density 
divided by (1 — m/2)(l — m^/4) in the case of the 
CI action. 

The considerations of this section wer based on 
the assumption of exact chiral symmetry. When 
using these results to interpret our data, we as- 
sumed that the chiral symmetry breaking effects 
in the CI action were small enough to be ne- 
glected. However, this assumption needs to be 
checked explicitly. 

Since the source is typically smeared, we mea- 
sured propagators from the smeared source to 
point sinks and smeared sinks. Appropriate ra- 
tios then allow us to recover e.g. unsmeared- 
unsmeared type propagators. 

^In our convention the experimental value of fr: is 
131MeV. 



Fig. ^ shows /tt as a function of the bare mass 
for Dci comparing a linear and a quadratic fit. 
For the extrapolation to the chiral limit we use 
the results from the quadratic fit. The chiral ex- 
trapolations of a/^ in Figs. 10 and 11 show ex- 
cellent scaling behaviour. 

4. Light hadron spectroscopy 

The parameters of our simulations (Table |l|) 
allow us to make a scaling study with both ac- 
tions in a fixed volume L « 1.3 fm at lattice con- 
stants a ~ 0.08, 0.10 and 0.15 fm and a finite vol- 
ume analysis at fixed a « 0.15 fm in boxes with 
L — 1.3 — 2.6 fm. We considered ~ 10 different 
quark masses and inverted the Dirac matrix with 
a multimass solver. Note that the hadron masses 
at different quark masses with other parameters 
fixed are strongly correlated. 

With this setup we could study different as- 
pects of light hadron spectroscopy close to the 
chiral limit. The subject is old, nevertheless we 
found surprises, results which we do not quite un- 
derstand and which require further study. 
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Figure 11. The (quadratic) extrapolations of /^r 
to the chiral hmit for different lattice spacings in 
dimensionless units (ro denotes the Sommer pa- 
rameter) exhibit simple scaling behaviour (results 
for Dci)- The dashed line corresponds to the ex- 
perimental (i.e. unquenched) value 131 MeV. 



4.1. Hadron spectroscopy in terms of 
hadronic scales 

Since our main interest lies in understanding 
the behaviour of the Dirac operators, it is natural 
to fix the scale within the hadronic sector itself. 
In Sect. 4.3 we shall connect these results to the 
typical gauge sector scale ro also. 

An obvious candidate to carry the scale could 
be the vector meson mass my at the quark mass 
where mps/'mv is equal to the experimental 
value. Identifying this vector meson mass with 
■nip « 770 MeV gives the cutoff and all the masses 
in MeV. This choice for the scale has the disad- 
vantage that it refers to the rho mass deep in 
the chiral limit and is plagued by relatively large 
statistical errors. For many purposes (like scal- 
ing and finite volume studies) it is better to fix 
the scale by the vector meson mass my(a;), where 
mps/mv = X fixed and chosen in a range where 
mvix) can be precisely determined. In our ana- 
lysis we have taken x — 0.75. A similar method 
was used earlier in the works |Q . 

Fig. ^ gives the vector meson mass mv as 
a function of the pseudoscalar mass mps in 
TOy(0.75) units. All our simulation results from 



Table are collected in this figure|. Indepen- 
dently of the lattice scale a, the size of the box 
(L = 1.3 — 2.6 fm) and for both Dirac operators 
the points lie on a universal curve. This observa- 
tion goes beyond our work: collecting data from 
other large scale simulations performed with dif- 
ferent Dirac operators (Wilson, clover, staggered 
fermions with or without improvement) the re- 
sults lie on this curve The point (0.75,1) of 
this curve is fixed by definition; the other points 
seem to be then independent of the cut-off in the 
range of recent simulations and in volumes down 
to L « 1.3 fm. 

Fig. |l^ shows the nucleon mass as a function of 
mps in my(0.75) units for the FP Dirac operator. 
The points forming the upper curve all refer to 
L w 1.3 fm, they represent a scaling test for a ~ 
(0.08 — 0.16) fm. Within the statistical errors 
no scaling violation is seen in this figure. The 
conclusion is the same for the CI operator. The 
large volume {L « 2.5 fm) results form the lower 
curve. These points were obtained on a coarse 
lattice with a 0.16 fm. 

For heavy quark masses we expect that the nu- 
cleon finite size effects will become smaller. We 
see indeed in Fig. |l^ that the two curves join for 
heavy quarks. This has an additional message. 
Since the upper curve scales, we might assume it 
is close to the continuum limit. This suggests that 
the large volume curve (for which we do not have 
a scaling test) is also, at least for heavy quarks, 
close to the continuum limit. 

Fig. illustrates that the scaling of the Dpp 
data in Fig. ^ is not 'built in' in the quantity 
studied. In this figure our large volume results are 
compared with the CP-PACS data from their re- 
cently completed analysis ||2^ . These results were 
obtained with the Wilson Dirac operator at four 
different lattice spacings (a = 0.05 — 0.10 fm). For 
the clarity of the figure we plotted the a = 0.05 fm 
and a = 0.10 fm data only. Even the a = 0.05 fm 
points are rather far from the continuum as esti- 
mated in [pl| . The difficulty of such a continuum 
extrapolation is illustrated in Fig. 



^Here and in the following, for the clarity of the figure, we 
suppress points with very large errors. 
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Figure 12. All the results on the vector and pseudoscalar masses of our simulations with different lattice 
spacings, volumes and actions. The scale is carried by the vector meson mass where mps/mv = 0.75. 
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Figure 13. The nucleon mass as a function of the pseudoscalar mass. The upper curve is a scaling test 
in a fixed volume. 
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Figure 14. Similar to Fig. 13, The large volume FP and CI results obtained on coarse configurations are 
compared with Wilson data and their continuum limit extrapolation [Ell . 
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APE plot where the FP and CI results are compared with Wilson |pi[ and improved staggered 
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Figure 16. The continuum extrapolation of the 
CP-PACS Wilson data H for 4 different x- 
values. The points on the far right are the FP 
a « 0.16 fm results. 



The four sets of points in Fig. |l^ refer to a; = 
iTips/mv — 0.75, 0.7, 0.6 and 0.4. For each set 
a linear extrapolation is done using the results 
at 4 different lattice spacings. The data and the 
extrapolated points at a = are taken from [2l| . 
The 4 points on the r.h.s. of Fig. 16 are the FP 
numbers at a « 0.16 fm for the 4 different a;- 
values. Further results are needed to make a firm 



conclusion from Fig. 14, 

Fig. |l^ is a standard APE plot where the large 
volume, coarse lattice FP a « 0.16 fm and CI 
a ^ 0.15 fm data are compared with the CP- 
PACS results and with the improved stag- 
gered fermion data obtained at a = 0.13 fm ||2q| 
and a = 0.09 fm ||2^ (2 points only). The coarse 
FP and CI data are significantly closer to the CP- 
PACS continuum extrapolated curve than the fine 
lattice Wilson or the improved staggered points. 
On the other hand, as we remarked above, the 
long continuum extrapolation of the Wilson data 
pil is not an easy task. 



4.2. Dispersion relation and the speed of 
light 

Another quantity which can be studied within 
the hadronic sector is the energy-momentum dis- 
persion relation. As earlier investigations show, 
E = E{p) is a cut-off sensitive quantity and the 
deviations from the continuum form (E'^(p) — 
iTL^)/p^ = can be large. Here c, the speed of 
light, should be independent of p and equal to 1 in 
our convention. For the FP Dirac operator at a « 
0.16 fm in a box with L « 2.5 fm the energy E{p) 
could be determined with relatively small errors 
in the momentum range \p \ = (0.5 — 1.0) GeV and 
the corresponding dispersion relation can not be 
distinguished from its continuum form 27 



Fig. shows the speed of light as a function 
of the quark mass for the vector and pseudoscalar 
meson as measured by the FP and the 3 over- 
lap steps augmented FP Dirac operators. For 
both hadrons the FP operator results show no 
deviation from 1, while the data with the over- 
lap augmented operator lie above 1 beyond the 
statistical errors. We find this result difficult to 
understand. Starting with the parameterized FP 
operator which approximates quite well a GW so- 
lution, we expected to produce small changes only 
by performing 3 overlap steps^. 

4.3. Connecting the hadronic and gauge 
scales 

In order to connect the hadron sector with the 
gauge sector it is sufficient to relate the reference 
scale TOy (0.75) we used above with a gauge theory 
scale, like tq The Sommer parameter has 

been measured for both actions close to the f3- 
range of our simulations [[lO|Jl4| . 

Fig. |l8| shows ror7iy(0.75) as a function of the 
lattice resolution. The 3 FP data on the coarsest 
lattice with different volumes L = 1.3 — 2.5 fm 
completely coincide showing that my (0.75) has 
no finite size effect within the small errors in these 
boxes. This is consistent with Fig. |lj. The 3 CI 
data on the coarsest lattice are consistent with 



^It is true, however, that in the spectrum of the overlap 
operator with the parametrized FP kernel even the low 
lying eigenvalues have a small shift along the circle relative 
to the spectrum of the kernel, see Fig(5.4) in the second 
reference in B. 
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Figure 17. The speed of light as obtained from the 
vector and pseudoscalar channels as a function of 
the quark mass. 



the FP number, but have larger errors. 

Relative to this coarsest point the a « 0.10 and 
0.08 fm data lie higher and indicate cut-off effects 
in roTOy (0.75). We have to remember though 
that measuring ro on a coarse lattice is a difficult, 
perhaps not even quite well defined problem. 

As Fig. |l^ shows, the prediction of the overlap 
augmented FP Dirac operator on the coarsest lat- 
tice is definitely inconsistent with that of the FP 
result. Since the Dirac operators FP and ov3 were 
used on the same gauge configurations (therefore 
they see the same a and Vq in a quenched the- 
ory) we have to conclude that one or both Dirac 
operators have clearly identified cut-off effects in 
mv(0.75) at a « 0.16 fm. For comparison we also 
included the Wilson data [|l) in Fig|l|. 

5. Pion scattering length 

The pion scattering length is an important ob- 
servable characterizing dynamical effects of the 
strong interaction and its ab initio calculation 
on the lattice is an important nonperturbative 
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Figure 18. ScaHng test for romy(0.75). 

test of QCD. In the literature several attempts 
to compute the scattering length with Wilson 
fermions | p^ , p9| , ^0pl] , |3^ and staggered fermions 
| |2^ , p9t can be found. However, calculations with 
a chiral Dirac operator which allows for a better 
control of the small mass region are still missing. 
In full QCD the scattering length is a quantity 
which vanishes in the chiral limit, while it is power 
divergent in the quenched theory and eventu- 
ally a study in the full theory is desirable. This is 
the initial stage of our investigation which we re- 
gard as a preparation for a simulation in the full 
theory when unquenched configurations become 
available. In this preliminary study we reuse the 
FP propagators computed for the spectroscopy to 
calculate the tttt S-wave scattering length in the 
1 = 2 channel. Although using these propaga- 
tors with only a single source limits the quality 
of the overlap with the pion scattering state the 
results are quite encouraging and a more refined 
investigation is in preparation. 

5.1. Method 

In order to calculate the Tr~^TT~^ scattering 
length from a Euclidean lattice simulation, we 
use Liischer's relation Q that relates the energy 
shift of the two pion state, AE{L), in a finite vol- 
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ume {V = L^) to the scattering length, oq, in the 
infinite volume limit to order 



AE = - 




where ci — —8.9136 and C2 — 62.9205 are nu- 
merical constants computed in |Q. Propagators 
generated for the spectroscopy study were used to 
construct the two pion state and to measure its 
energy, E = 2 + Ai?. In other words, a single 
Gaussian source was used to create the two pions 
on the gauge fixed configurations. We also ap- 
plied the same Gaussian smearing at the sinks in 
this study. We worked on the 8^ x 24 and 12^ x 24 
lattices at a = 0.16 fm. All the scattering lengths 
shown in Fig. |l] are for i > 4 except for the 
two lightest masses for L = 12 a where ttIt^ L is 
only larger than 3. 

5.2. Analysis 

We extract the energy shift of the two pion 
state in a finite volume by considering ratios of 
two pion correlation functions and single pion cor- 
relation functions. We follow the procedure and 
notation of Ref. |^ and construct the "direct" 
(D) and "crossed" (C) contractions of 

l^Oi{xi,t)Y,02{x2,t)S{t = Q)S{t = Q)\ , 

(11) 

where O annihilates a 7r+ and S are the gaussian 
sources for the pions. Note that D and C actually 
refer to the ratio after dividing by two single pion 
correlators. In practice we fit the numerator and 
denominator separately. The different diagrams 
(contractions) of ratios are labeled, 

D{t) ~ glue exchange between the two pions 
C{t) ~ quark exchange between the two pions 

where the combination of interest (/ = 2) is 
D{t) — C{t) whose large t behaviour is approx- 
imately Ze^^^*. We have also expanded the 
exponential to leading order in AEt in the fit, 
but there were no differences within the errors. 

Periodic boundary conditions were used in the 
time direction so that we must take into account 
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Figure 19. The numerator of the D — C function. 

a backwards propagating pion. The fitting form 
of the two-particle correlator (the numerator of 
D{t) - C{t)) is 

Z^^(e"(^™''+^-^)*-t-e"(^™''+'^-^)('^"*)+ze"'"''^) . 

(12) 

We simultaneously fit the pion propagator enter- 
ing the denominator for Z^r and m^, 

Z^(e-™-* + e-"-(^-*)) , (13) 

resulting in the following five fitting parameters: 
TOtt, ^tt, AE, Ztttt and z. We perform a fully cor- 
related fit to the two correlation functions. In the 
two figures that follow, we show the fits for the 
following set of parameters, 

Pfp = 3.00, a « 0.16 fm 
Ns ^ 12, Nt = T = 24 
mir/mp ~ 0.41, 200 configurations 

which is representative of all the other fits. 

In Fig. H^, we show how the extracted ao de- 
pends on the initial time slice included in the fit. 
We see that nearly any time slice from t = 2a 
gives the same results (with reasonable val- 
ues) within the errors. This may be somewhat 
surprising considering that our two pion sources 
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The dashed Figure 21. Summary of our results and other re- 

cent calculations. 



are on top of each other. We take the value for the 
fit starting at tmin = 3 a in this particular case. 
There were substantial difficulties in getting re- 
sults for the other lattice spacings as the tempo- 
ral extent of the box was not as long and/or short 
of statistics. In the following, we mainly discuss 
the results from the lattice described above. 

5.3. Results and discussion 

Using the parameterized fixed-point action, we 
have extracted the 7r+7r+ S-wave scattering length 
at threshold with the propagators generated for 
spectroscopy studies. We were able to obtain re- 
sults for the l3pp = 3.0, 12^ x 24 lattice and for 
some of the larger masses on the 8'^ x 24 lat- 
tice. In Fig. we plot the dimensionless quan- 
tity (as was used in Ref. |3l| ), ni^aQ/niT^, against 
rriTr/mp. We also include in this plot some recent 
results from other groups as well as quenched chi- 
ral perturbation theory |p3| and chiral perturba- 
tion theory predictions ||35| . Our raw results on 
a coarse lattice are in reasonable agreement with 
Wilson/clover action results on finer lattices or 
continuum limit extrapolated numbers. Since we 



mostly only have results from a single coupling, 
no real scaling tests could be performed for qq. 

However, agreement between two of the coarser 
lattices in limited cases for rriTrL > 4 and mT^/mp 
ratios larger than 0.6 indicate that they are small 
(also found in spectroscopy studies). The difficul- 
ties of extracting the numbers from the other lat- 
tices are likely due to contamination from higher 
momentum pion states. These effects can be re- 
duced by introducing another source, which is 
currently under study. As remarked above some 
zero mode effects were observed in the pseu- 
doscalar propagator in the spectroscopy study for 
fnps/'fnv mass ratios less than 0.5. These finite 
size effects as well as other systematic error stud- 
ies are work in progress. 

6. Summary and conclusions 

We have reported here several calculations of 
quantities relevant for light quark phenomenology 
using two different Dirac operators which are ap- 
proximate solutions of the Ginsparg- Wilson equa- 
tion. 
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Both actions, beyond having good chiral be- 
haviour, exhibited significantly smaller cut-off ef- 
fects at a ^ 0.15 fm than the standard Wilson 
action at a '--^ 0.05 fm and the improved-staggered 
action at a ~ 0.09 fm on hadronic quantities that 
we have surveyed. 

We believe that the improvement more than 
compensates the apparent increase in the cost of 
simulating our operators and that a study of dy- 
namical simulations with these actions is a rele- 
vant challenge for the future. 
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